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The  problem  of  unsteady  motions  of  a thin  airfoil  under  the  free  surface 
of  a weightless  fluid  of  finite  depth  is  discussed  in  the  linear  formulation. 
The  analogy  to  flow  around  a cascade  with  an  infinite  number  of  special-shape 
airfoils  is  employed.  The  method  of  discrete  vortices  is  used  to  solve  the 
latter  problem. 


Let  us  consider  small  harmonic  oscillations  of  an  infinitely  thin  airfoil 
of  length  2 c in  a plane-parallel  flow  of  a weightless  incompressible  fluid  with 
constant  undisturbed  velocity  Vq,  bounded  by  a free  surface  above  and  by  a solid 
bottom  below. 


The  origin  of  the  body  axis  system  will  be  taken  at  the  center  of  some  mean 
position  of  the  oscillating  airfoil.  The  Ox  axis  will  be  directed  against  the 
incident  flow,  and  the  Oy  axis,  vertically  upward.  The  mean  submersion  depth 
of  the  airfoil  under  the  free  surface  will  be  denoted  by  h,  and  the  mean  distance 
from  the  bottom,  by  hg.  In  the  case  of  harmonic  oscillations  of  frequency  p, 
the  disturbance  potential  may  be  represented  as 

<t>(x,  y,  t)  = 4>(x,  y)e±pt. 

Here  and  below,  only  the  real  part  of  the  expressions  containing  the  time  factor 
in  exponential  form  should  be  used. 


For  the  reduced  potential  <J>(x,  y) , a linear  boundary  value  problem  will  be 
obtained  on  the  basis  of  the  accepted  assumptions. 


It  is  necessary  to  find  the  solution  of  Laplace's  equation 


satisfying  the  boundary  conditions 
4>x  = 0 for  y = h (2) ; 


In  addition,  according  to  the  requirement  of  the  Zhukovskiy-Chaplygin 
hypothesis,  the  velocity  boundedness  condition  in  the  trailing  edge  of  the 
airfoil  should  be  fulfilled  for  y = 0 and  x = -c.  , 


♦Numbers  in  the  right  margin  indicate  pagination  in  the  original  text 
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The  boundary  value  problem  (l)-(5)  will  be  solved  by  the  method  of  singu- 
larities, using  the  layer  of  bound  vortices  along  a segment  of  the  Ox  axis 
(y  **  0;  -c  < x 4 c) . The  strength  of  the  free  vortices  of  the  airfoil  and 
wake  will  be  determined  on  the  basis  of  the  Birnbaum-Kussner  unsteady  flow 
pattarn.2 

It  is  easy  to  show  that  the  problem  uqder  consideration  is  reduced  to  the 
problem  of  unsteady  flow  past  some  special-type  infinite  lattice.  Such  an 
analogy  is  based  on  the  fact  that  condition  (2)  on  a free  surface  in  vortex 
flow  can  be  satisfied  by  placing  a vortex  with  a strength  of  the  same  magnitude 
and  sign  at  a distance  h above  the  free  surface. 


To  fulfill  the  condition  on  the  solid  boundary  (3) , a 
vortex  with  a strength  of  equal  magnitude  but  opposite 
sign  must  be  positioned  symmetrically  relative  to  this 
boundary.  Thus,  for  a single  bound  vortex  in  a flow  of 
fluid. of  finite  depth,  we  arrive  at  the  vortex  chain 
shown  in  Fig.  1. 

We  will  determine  the  vertical  velocity,  caused  by 
the  vortex  chain  in  question,  at  some  point  of  the  flow 
M(x,  y) . If  the  vortex  strength  is  constant, 
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At  the  points  of  the  Ox  axis,  we  have 
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Using  the  equality 
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we  sum  up  expression  (6)  and  obtain 
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For  a variable  vortex  strength  T = T(t),  the  velocities  in  the  plane  of 
flow  are  induced  not  only  by  point  vortices,  but  also  by  their  vortex  sheets 
with  a linear  vortex  strength 

ir^ntn 


Then  on  the  basis  of  the  Birnbaum-Kiissner  hypothesis 
given  by 

v,(x)  ~ (v™  + 

where  0 r 

= vM (x)  — k f sinfo'oro(j;  — s^df; 


, the  vertical  velocity  is 


is  the  Strouhal  number 


To  find  the  nonsteady  aerodynamic  characteristics  of  such  an  airfoil,  it 
is  necessary  to  determine  the  strength  of  the  bound  vortices  of  the  airfoil.  On 
the  basis  of  the  no-flow  condition  (4)  of  a continuous  vortex  layer,  expressions 
(9)  and  (10)  make  it  possible  to  write  an  integral  equation  for  finding  the  unknown 
strength. 


The  nonsteady  aerodynamic  characteristics  of  the  airfoil  are  conveniently 
expressed  in  terms  of  the  coefficients  of  rotary  derivatives.^- 

The  unknown  vortex  strength  of  a nondeformable  airfoil  will  be  represented 
in  the  linear  approximation  as  follows: 


where  a is  the  angle  of  attack;  w is  the  absolute  angular  velocity  of  the  body 

axis  system,  a = ; w = the  coefficients  Y^2)  are  referred  to  ai 

the  rotary  derivatives  of  the  vortex  strength. 

The  kinematic  characteristics  of  the  airfoil  will  be  assumed  to  change  in 
accordance  with  the  laws 


In  view  of  expressions  (12) , the  vortex  strength  (11)  will  be  written  as  follows 


Y =>  yWv"  + 


We  transform  the  no-flow  conditions 

vy  = -VQa  +.ojx 

by  means  of  formulas  (12)  to  the  form 
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To  set  up  the  integral  equation,  it  is  necessary  to  calculate  the  vertical 
velocity  caused  by  the  adjacent  surface  layer  with  strength  (13),  and  subordinate 
it  to  the  condition  (14).  The  equation  thus  obtained  turns  out  to  be  singular 
and  in  addition,  to  contain  a logarithmic  type  singularity  - (cikx) . 

Separating  the  real  from  the  imaginary  parts  in  the  integral  equation,  and 
in  the  imaginary  parts,  the  terms  containing  equal  coefficients  or  we 
obtain  two  systems  of^ integral  equations,  which  we  write  in  dimensionless  form 

j Yi0  (s)  ^ ( x — s)ds  — k j Yl2,(^2)  (x  — J)  ds  = — 2n; 
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* = T;  S = T;  *-£•. 

The  method  of  discrete  vortices  (1)  will  be  used  to  solve  the  system  of 
integral  Equations  (15)  and  (16) . 

Mathematically,  this  method  signifies  the  replacement  of  definite  integrals 
within  the  limits  (-1,  +1)  by  formulas  of  mechanical  quadratures  with  nodes  at 
4i  - 1 

points  s.  *=  -1  + J , — and  an  approximate  fulfillment  of  the  equations  of  the 
J 4i  - 3 

system  at  some  points  Xi  _1  + 4N  close  to  nodes  s ^ . Such  a shift  of 
variables  in  singular  integrals  of  the  form 

J y{s)---~fds,  — 1<.T<1 

makes  it  possible  to  avoid  kernel  singularities  when  replacing  an  integral 
equation  by  a system  of  algebraic  ones. 

In  the  indicated  arrangement  of  points  x^  and  sj , a boundedness  of  the 
function  Y(x)  is  reached  when  x ->  -1;  when  x -*■  1,  it  has  a singularity  of 
integrable  order.  For  other  classes  of  functions  Y(x) * other  layouts  of 
points  x^  and  sj  can  be  proposed. 
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The  aerodynamic  coefficients  will  be  represented  in  terms  of  the  total 
rotary  derivatives 
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which,  according  to  the  formulas 
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C„--=  f Y (*)dx;  C„  — [xy(x)dx. 


are  expressed  in  terms  of  the  distributed  rotary  derivatives  y^1^ » • • • , y(2) 


When  the  method  of  discrete  vortices  is  used,  the  integrals  will  be  re- 
placed by  the  finite  sums 
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where  N is  the  number  of  bound  vortices. 


u T£e  valldlty  criterion  of  the  calculations  is  the  requirement  CW  = Ca  and 
Cy  * Cm*  which  is  a corollary  of  the  reversibility  theorem.1 
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FIGURE  2 


FIGURE  3 
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Figures  2 and  3 are  plots  of  the  change  in  the  functions  ipa  = — ^0  and 


f', 


Cmhhg 


a 


m 


ca 

Cn|co 


in  steady  motion  (k  = 0)  with  the  relative  distance  of  the  thin 


•airfoil  from  the  bottom.  It  follows  from  these  figures  that  as  the  distance 
increases,  the  aerodynamic  coefficients  of  the  airfoil  approach  the  values  of 
the  corresponding  characteristics  during  motion  under  the  free  surface  of  a 
fluid  of  infinite  depth.  As  the  distance  of  the  airfoil  from  the  bottom  de- 
creases, its  aerodynamic  characteristics  increase  and  approach  the  corresponding 
values  for  motion  near  a solid  boundary. 
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FIGURE  4 


Figure  4 shows  the  change  in  the  position  of  the  pressure  center  of  a 
thin  airfoil,  a change  barely  noticeable  near  the  free  surface.  As  the  solid 
boundary  is  approached,  the  pressure  center  shifts  appreciably  toward  the  mid- 
section  of  the  airfoil. 


The  numerical  results  obtained  for  steady  flow  at  relative  submersions  and 
distances  h >0.25,  hg  > 0.25  are  in  good  agreement  with  the  results  of  A.  N. 
Panchenkov's  calculations,^  obtained  by  using  the  asymptotic  small-parameter 
method.  At  small  submersions  and  distances,  the  method  described  makes  it 
possible  to  obtain  results  corresponding  more  accurately  to  the  initial  equation. 
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